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We revisit the problem of the electromagnetic Green function for homogeneous hyperbolic media, 
where longitudinal and transverse components of the dielectric permittivity tensor have different 
signs. We analyze the dipole emission patterns for both dipole orientations with respect to the 
symmetry axis and for different signs of dielectric constants, and show that the emission pattern 
is highly anisotropic and has a characteristic cross-like shape: the waves are propagating within a 
certain cone and are evanescent outside this cone. We demonstrate the coexistence of the cone- 
like pattern due to emission of the extraordinary TM-polarized waves and elliptical pattern due to 
emission of ordinary TE-polarized waves. We find a singular complex term in the Green function, 
proportional to the function and governing the photonic density of states and Purcell effect in 
hyperbolic media. 
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I. INTRODUCTION 

Hyperbolic medium is a particular case of uniaxial 
anisotropic dielectric medium where the main values of 
the permittivity tensor e have opposite signs [1], and the 
isofrequency surface of extraordinary waves is a hyper- 
boloid. The recent attention of researchers to the hy- 
perbolic metamaterials stems from their unique optical 
properties allowing negative refraction, hyperlensing and 
cloaking phenomena [2|-l4( . They are also very promising 
for quantum nanophotonics Q, because the density 
of states, determined by the hyperboloid area, is diver- 
gent. This means an infinite spontaneous decay rate of a 
quantum emitter embedded in a hyperbolic medium, i.e. 
infinite Purcell effect [7[. In the realistic case, the radi- 
ation rate is limited by certain cutoffs in the wavevector 
space j^llj, however, experimental observation of radia- 
tive enhancement is still possible [l2| . 

A general hyperbolic medium can be characterized by 
a dielectric permittivity tensor 




Ree± Res\\ < . 



(i) 



Resulting isofrequency surface of extraordinary waves is 
a hyperboloid, 



o 2 - 



(2) 



Two different types of hyperboloids are possible for ey < 
0, £_l > (see Fig. [TJi) and for s\\ > 0, e± < (see 
Fig. Eh). 

Hyperbolic medium can be realized in several ways. 
First realization has been reported for magnetized 
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FIG. 1: Schematic illustration of the isofrequency surfaces in 
the wavevectors space for the hyperbolic medium with su < 0, 
e± > (a) and with s\\ < 0, e± > (b). 



plasma in microwave spectral range [l3|- Under a very 
strong static magnetic field applied along z axis the 
plasma is described by the permittivity with [l3|, QjJ 

uul 

£± = 1, £|| = 1 o • 

Thus, the frequencies uj below the plasma frequency uj p 
correspond to the hyperbolic regime with e\\ < 0. 

The interest to the hyperbolic medium is now revived 
and rapidly increases due to its successful realization us- 
ing artificial photonic structures, metamaterials 0, flBI - 
Il8j . In particular, the layered structure composed of al- 
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ternating dielectric and metallic slabs is described by ef- is given by 
fective permittivities 



di + d 2 



e 1 x di 



~ X d 2 



d 2 



where e±, £2, d 2 are dielectric constants and thick- 
nesses of slabs, respectively. Since dielectric constants of 
the metals in optical frequency range are negative, it is 
possible to adjust thicknesses of slabs in order to obtain 
either e± < or en < [19|, [20|. Realizations of hyper- 
bolic regime have been also rep orted for metamaterials 
based on nanorod arrays [l6|,[2l|,[22| and for graphite [23| . 

The ongoing studies of hyperbolic medium raise the 
demand for the general theoretical formalism. The most 
rigorous description of the optical properties of arbitrary 
photonic structure is provided by its tensor Green func- 
tion [24[, determined from 



(VxVxI- k 2 e)G(r) = 47rk 2 I5(r) , 



(3) 



G(q) = 47r/c 2 < e±e\\i ■ 



q (g> q 



k 2 J q 2 s\\ +q 2 L e± - k 2 e±£\\ 



(q x z) (g) (q x z) 



g 2 — /c 2 £x 



g 2 £|| - fc 2 e_L£|| 



(5) 



Here q\\ = ^ and = ^ + The symbol (g) de- 
notes the direct product, [a 6]^ = dab p. The poles in 
Eq. (|5|) determine the dispersion equations of the electro- 
magnetic waves. Two poles correspond to extraordinary 
(TM) waves, with magnetic field H perpendicular to z 
axis, and to ordinary (TE) waves, with E • z = and 
q = k^Jel- 

Both direct and reciprocal space methods give the 
same following result for the Green function: 



where / is the unit tensor and k = lj/c. Green function 
for medium with uniaxial permittivity tensor is presented 
in a number of works [l],[25-33]. Most general form of the 
result is given by Chen in Ref. [26| in dyadic form and by 
Savchenko in Ref. [32[ in Cartesian form. However, de- 
spite the impressive amount of the researchers done, we 
are not aware of any comprehensive study of the Green 
function in hyperbolic case. Moreover, all existing works 
except [29|, |31( neglect the singular term in the Green 
function, which, as will be shown below, is crucial for the 
description of the photonic density of states and Purcell 
enhancement in hyperbolic regime. Here we present a 
general theory of the Green function and dipole emission 
in hyperbolic medium. We analyze both types of hyper- 
bolic medium, illustrated on Fig. [TJ and both axial and 
transverse dipole orientations. 

The rest of the paper is organized as follows. Sec. [XT] 
outlines the Green function calculation. Singular term 
in the Green function is discussed in Sec. IIIII Sec. IIVI 
presents the calculated emission patterns of point dipole. 
Main paper results are summarized in Sec. [Vj Auxiliary 
expressions are presented in Appendices. 



II. GREEN FUNCTION CALCULATION 

Green function (j3j) may be calculated either in real 
space via the methods of operators [30| or in wavevector 
space via Fourier analysis [26]. The Fourier representa- 
tion of G(r), defined from 



d 3 q 
(2tt) 3 



iqr 



G(q) 



(4) 



G(r) 



a ikr e 



e + V<8> V 



ikr Q 



^ikv ^ikv 

(r x z) 2 



n ikr e 



k 2 e ± 



ik 



(r x z) (g) (r x z) 



i — z & z 



(r x z) 2 
2(r x z) ® (r x z) 



(6) 



(r x z) 
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where 



e = e±e\\i 1 = e\\ (I - z <g> z) + e±z ® z, 
r e = \/r(er) = J £\\(x 2 + y 2 ) +£±z 2 , r Q = 



e±r. 



This result was obtained without any assumptions for 
signs of real of parts e± and e\\ and is applicable to case 
of hyperbolic medium. Imaginary parts of e± and su and 
of the square roots in expressions for r e and r Q should be 
positive (we assume the time dependence e~ luJt ). 

Eq. ([6|) is not the final result for the Green function 
yet, because the derivative in the term V Ve lkre /r e is 
not evaluated. Some authors [30,[33[ leave the expression 
for G(r) in the form (|6]) without taking this derivatives. 
However, this term produces singularity in G(r). The 
nature of this singularity is the same for the well-known 
identity 
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-47r5(r) . 



Calculating the derivative we obtain the final expres- 
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sion of singular Green function: 
1 



ikr e 



G(r) = G sing (r) 
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ikr e 



k 2 7 



n ikr e 



n ikr Q 



ikr P 



ikr a 



n ikr e 



(er) <g> (er) 
(r x z) (r x z) 
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n ikr e 
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(r x z) 2 
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where G S ing(r) stands for the singular contribution 



G S ing(r) 



(V<8) V) 



-47T 



sing 



(2tt) 3 



2qr _ 



q q 



(8) 



In general case we failed to obtain a closed answer for 
Gging via ordinary functions and (5-function. Eq. (|8j) 
should be understood instead only as generalized func- 
tion, i.e., only its convolutions with ordinary test func- 
tions are relevant [34|. Singular Green function is es- 
sential in the hyperbolic case since it solely accounts for 
the diverging Pur cell factor. Eqs. (O,® are the central 
result of this work. It is valid for both hyperbolic and el- 
liptic media with arbitrary signs of real parts of dielectric 
constants. 



III. SINGULAR TERM OF GREEN FUNCTION 

Singular term in the Green function should be treated 
with extra care [35-37]. Functional (j8j) may be explicitly 
written out in only isotropic medium with s\\ = e± = 



Gsing(r) 



47r r ( 



-S(r) . 



(9) 



In general anisotropic case we have not found a closed 
form answer. Previous expression for this functional in 
anisotropic medium, obtain by Weiglhofer [28|, [29| , 



(10) 



is obviously wrong since it has been obtained as a naive 
generalization of the result [39] 



47T 

3e 



S(t) 



(11) 



for the isotropic medium. Despite the fact that func- 
tional ([TT]) is frequently used [39| , it gives different results 
than Eq. (j9]) when applied to functions, non-analytical in 
the point r = 0, such as r (g) r/r 2 . This discrepancy is 



discussed in Ref. [38] in details. Since angular averaging 
of r (g>r equals to 1/3, both functionals ([TT]) and ([9]) may 
give the same result for some test functions. However, 
straightforward generalization of Eq. ([TT]) to Eq. ([TP]) is 
not valid since not all rules of normal calculus may be 
applied to generalized functions [34j . 

Explicit form of the functional (|5} in anisotropic case 
can be obtained within the space of the test functions, 
analytical in the point r = 0: 



^sing(r) 



e\\ - e± 



+ arctan(w£ii/£_L — 1) 



\V £ \\/ £ ^ - 1 
(x®x + y(g)y) 



z 09 z 



(12) 



The equivalence between Eq. ([T2]) and Eq. ([8]) for test 
functions, analytical in the point r = is directly shown 
via real space integration in spherical coordinates. Let 
us demonstrate it for zz component: 



J d 3 r[G sing ] zz (r) = j d 3 



d 2 



r<R 



r<R 



dz 2 yfeZr e 



dS ■ z- d 1 



r=R 



dz ^Jelr e 



sin v cos 



(en sin 2 6 + e ± cos 2 Of/ 2 



(13) 



Performing the integral over 6 we obtain expression, 
equal to J d 3 r[G S i ng ] zz (r) which finalizes the proof. 
Eq. ([T2]) also follows from the wavevect or- space represen- 
tation in Eq. ©. In the limit e± = £\\ Eq- ([T2]) reduces to 
Eq. ([TT]) . Obviously, functional ([T2]) is far more complex 
than (wrong) Eq. ([TO]) . Still, it is valid for narrower set 
of test functions than Eqs. (|8]),(|9]). 

One may think that this singular terms are of purely 
mathematical interest. Counterintuitively, they may con- 
trol such observable quantities, as photonic density of 
states and Purcell factor. In particular, in the case of 
lossless hyperbolic medium (s\\s± < 0, linen, £j_ — » +0) 
the singular part of the Green function acquires non-zero 
imaginary part, 



Im G S ing(r) 



Nl) 3/2 

\e ' 



S(r) 



(x (g) x + y (g) y) 



(14) 



The singulary in Eq. ([T4]) is a direct consequence of the in- 
finite density of TM modes Eq. (|2]). Green function where 
this singularity is omitted has wrong analytical proper- 
ties and can not be used for nanophotonic applications 




Re e± > 
Re e, < 




t 



-1 1 



FIG. 2: Dipole field \ 3 E z (0,y, z) in hyperbolic medium with e± = 1, ey = —3 (a) and with s_l = 1 + 0.2z, ey = — 3 + 0.2i (with 
losses) (b). Dipole moment is parallel to the axis of anisotropy z. 



such as Purcell factor calculation. Generally, the Pur- 
cell factor of the embedded light emitter oriented along 
n direction can be found [40| as 



/=2fc3lmn.G(0)n. 



(15) 



Due to the singularity in the local density of states 
Eq. ([T5]) provides diverging Purcell factor. As has been 
indicated in our previous work [8[, this divergence is 
smeared out for finite size emitter, characterized with 
spatial distribution 3>(r), normalized as J d 3 r&(r) = 1. 
In case of the semiconductor quantum dot 3>(r) is pro- 
portional to the exciton envelope function [4l|. For dis- 
tributed source Eq. ([T5]) is replaced by 



/ = • J <PrcPr'$(r)$(r')G(r - r')h . (16) 

In the case of isotropic source distribution, 3>(r) = 3>(r), 
the Purcell factor is readily evaluated using (fTl| and is 
proportional to the cube of the ratio of the wavelength 
and the size of the source. 



IV. DIPOLE EMISSION PATTERN 

Green function (|3j) allows to find electric field for arbi- 
trary distribution of polarization P(r), 

E(r) = J d s r f G(r - r')P(r') . (17) 

For the point dipole p one has P(r) = p5(r) and the 
electric field is given by 



In the following subsections we consider two principal 
cases of orientation of dipole moment p with respect to 
the anisotropy axis z. 

It should be noted that the regular part of the dipole 
field may be easily obtained without knowledge of the 
dyadic Green function ©• Electromagnetic field in the 
uniaxial medium can be decomposed into two parts: TM- 
field, where H • z is zero, and TE-field, where E • z is 
zero [25| . Solutions for these fields in hyperbolic medium 
can be found separately via corresponding solutions in 
vacuum using of anisotropic scaling method introduced 
by Clemmow in Ref. [42|. The key point of this method 
consists in appropriate scaling of space as well as fields 
and polarizations from Maxwell's equations in vacuum to 
obtain expressions for fields and currents in anisotropic 
medium. Application of Clemmow's method for hyper- 
bolic medium gives the following scaling rules: 
(i) TE polarization 



E(r) = G(r) • p . 



(18) 



E(r) = e 1/2 E (e 1/2 r), 

H(r) = /I^Ho(e 1/2 r), 

P(r) = VeI^e 1/2 Po(e 1/2 r), 

(ii) TM polarization 

E(r) =E ( V / iIr), 
H(r) = v^IHo(V^Ir), 
P(r) = £±P (V^Ir), 

where Eq, Ho, Po are vacuum solutions and 

£ 1/2 = - z ® z) + y'Ejjz <g> z, 

e 1/2 = y/ejtfl - z <g) z) + y/slz ® z. 



(19) 
(20) 
(21) 

(22) 
(23) 
(24) 
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FIG. 3: Dipole field \ 3 E Z (0, y, z) in hyperbolic medium with e± = 
losses) (b). Dipole moment is parallel to the axis of anisotropy z. 



-1, £|| = 3 (a) and with e± = — 1 + 0.2z, e\\ — 3 + 0.2i (with 



Decomposition of the dipole polarization P into the 
TM/TE parts is described in [II S3- 



A. Dipole parallel to the axis of anisotropy 

Here we consider the case p = pz. Taking into account 
that 

i • z = s±z, (er) • z = e±r • z, 
(r x z) • z = 0, (I — z (g) z) • z = 0, 

we obtain from Eqs. (0, (fT8|) the following result for the 
regular part of the electric field: 



[E(r)] reg =^ 



p ikr e 
Akr e 



1 



1 



1 



ikr e 
3 



k 2 r 2 
3 



ikr e k 2 rl 



z (25) 
(r -z)(er) 



Cartesian representation of this expression is presented in 
lAl We see that the axial dipole emits only TM-polarized 
(extraordinary) waves. Calculated cross-section of the 
electric field in the yz plane is presented in Fig. [2j Fig. [3] 
for two both cases of hyperbolic medium, illustrated in 
Fig. Q] Reey < 0,Re£ ± > and Res\\ > 0,Ree ± < 0, 
respectively. Electric field pattern has a distinct cone-like 
shape: the waves are emitted only within the polar angles 
satisfying en sin 2 + e± cos 2 > 0. This radiation pat- 
tern, characteristic for hyperbolic medium, survives even 
when the losses are introduced, see Fig. [2Jd and Fig. [3b. 
An interesting feature of Eq. (jIVBjh revealed in Fig. [2j 
is that |^(0, 0, z)\ is not zero for axial dipole orientation 



and decays as l/\z\. This means that in the anisotropic 
medium the field is nonzero even along the direction of 
the dipole z. 



B. Dipole orthogonal to the axis of anisotropy 

Here we consider the case p = py. Taking into account 
that 

e-y = £||y, (er)-y = £||r-y, 
we obtain from (j7j), ([18]) the following result: 



ikr e 



[E(r)] reg =^ 2 



■pk' 
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-pk 2 l^/el 
ipk 



ikr e 

3 
ikr e 



k 2 r1 
3 



(r -y)(er) 



ikr a 



ikr e 



To 

ikr Q 



[(r x z) -y](r x z) 



£± r, J (r x z) 2 

2[(r x z) -y](r x z) 



y/e± (r x z) 2 



(r x z) 2 



(26) 



Cartesian representation of the last expression is pre- 
sented in[Bl For this geometry both TE and TM polarized 
waves are emitted. Equation (|26|) includes the denomi- 
nator (r x z) 2 = x 2 + y 2 which turns to zero at the line 
x = y = 0. However, careful analysis of Eq. ([26]) ensures, 
that electric field is continuous at this line since diverging 
TE and TM wave contributions cancel each other. Calcu- 
lated emission pattern is presented in Figs. I4)5l for differ- 
ent signs of dielectric constants. Most interesting results 



6 



(a) 5 




(b) 5 



2 - 



-3 - 



Re e± > 
Re e, < 




=> 



-5 -4 -3 -2 -1 1 2 3 4 5 



n 



U 



FIG. 4: Dipole field \ 3 E y (0, y, z) in hyperbolic medium with e± = 1, ey 
losses) (b). Dipole moment is orthogonal to the axis of anisotropy z. 
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FIG. 5: Dipole field \ 3 E y (0, y, z) in hyperbolic medium with s± = —1, s\\ = 3 (a) and with s± — —1 + 0.2i, ey = 3 + 0.2i (with 
losses) (b). Dipole moment is orthogonal to the axis of anisotropy z. 



are manifested for Re en < 0,Re£^ > (Fig. 3]), when 
the electric field is a distinct superposition of the conic 
patern due to the TM waves emission and elliptic pattern 
due to TE waves. In the second case Re en > 0, Re e±_ < 
(Fig. [5]) the TE waves contribution leads just to the spa- 
tial modulation of the conic radiation pattern. Similar 
to the case of the axial dipole orientation, the far field is 
present even along the dipole direction. 



Interesting results are also revealed in the Poynting 
vector distribution, shown in Figs. I6|71 The Poynting 
vector is found as S = c/(47r)ReE x H*, explicit ex- 
pressions for magnetic field H are presented in lA|B[ The 
Poynting vector pattern inherits the conical shape from 
the electric field. Its largest values are achieved at the 
conical surface r\ = 0. Interesting features are observed 
for the Poynting vector distribution in case of perpen- 
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FIG. 6: Poynting vector 4tt/ (cA 6 )S(0, y, z) in hyperbolic medium with e± — 1 + 0.22, s\\ — — 3 + 0. 2i (a) and with e± — —1 + 0.22, 
£\\ = 3 + 0.2i (b). Dipole moment is parallel to the axis of anisotropy z. 
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FIG. 7: Poynting vector 47r/(cA 6 )S(0, y, z) in hyperbolic medium with e± — l + 0.2z, s\\ — — 3 + 0.2i (a) and with e± — —1 + 0.22, 
£\\ = 3 + 0.2i (b). Dipole moment is orthogonal to the axis of anisotropy z. Inset in panel (a) illustrates the Poynting vector 
distribution in the xz plane for y = 30 A. 



dicular dipole orientation and Reen < 0, illustrated on 
Fig. [7^. In the region outside of the cone, r\ < 0, the 
Poynting vector in the yz plane is directed almost along 
the z direction. The resulting pattern looks as if there is 
a line source located at y = x = 0, which violates the en- 
ergy conservation condition. However, this first impres- 
sion is wrong. Close inspection of the Poynting vector 
distribution reveals that the line y = x = is a saddle 



point of the Poynting vector: the energy enters the line 
along x direction and leaves along z one. This behavior is 
illustrated in the xz cross-section of the Poynting vector 
distribution, shown in the inset of Fig. [7^. 



One can also prove this result analytically. Neglect- 
ing the evanescent terms oc exp(— k\r e \) we find following 
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expression for the Poynting vector for vanishing losses 



S = 



47rr 3 / 2 (x 2 + ^ 2 ) 2 v^T 
x ^k[k 2 e±x 3 (x 2 + y 2 ) + x(x 2 - y 2 )} 
+y[k 2 e±x 2 y(x 2 + y 2 ) + 2x 2 y] 
+z[k 2 e±x 2 z(x 2 + ?/ 2 ) + z(r 2 + 



(27) 



Direct differentiation of Eq. (|27|) demonstrates that V • 
S = 0, i.e., energy conservation law is fulfilled. In the 
limit \x\ : \z\ <C y, 1/k we find S ex (zz — xx)/|?/ 4 |, which 
confirms existence of the saddle point in the Poynting 
vector pattern. 



Appendix A: Cartesian representation of the field of 
dipole parallel to the anisotropy axis 



Cartesian components of the electric field Eq. (|IVB|) 
of the dipole oriented along the anisotropy axis z read 



E x -pk 2 e\\y/el- 



E y -pk 2 s\\^T[- 



E z = pk 2 £\\^/eZ- 



ikr e 
ikr e 



2 , 2 

x + y 



ikr e k 2 r 2 



ikr e k 2 r 2 
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yz, 



(Al) 



x 2 + y 2 _ *^± z 2 

ikr P ' k 2 r 2 ) \ ' s\\ 



V. CONCLUSIONS 



where r e = ^/e\\(x 2 + y 2 ) + s±z 2 . Similar result has 
been obtained in Ref. [32]. Magnetic field can be found 
as H = -{i/k)V x E: 



We have presented a general theory of the dipole 
emission in homogeneous hyperbolic media. Using both 
Fourier space approach and electromagnetic scaling, we 
have obtained a general expression for the electromag- 
netic Green function, and demonstrated that the emis- 
sion pattern is highly anisotropic. For dipole orienta- 
tion parallel to the symmetry axis, only TM-polarized 
waves are excited and the emission pattern has a cone- 
like shape with propagating waves present only within 
the cone en sin 2 + e± cos 2 6 > 0, where 6 is the po- 
lar angle. In case of the perpendicular orientation, the 
electric field is given by a sum of TE-polarized and TM- 
polarized contributions, so the waves can propagate also 
outside the cone. We have revealed the importance of 
a singular term in the Green function, and have demon- 
strated that it is of crucial importance for calculation 
of the radiative rate of light source embedded in a hy- 
perbolic medium. In the conventional case, the singular 
term proportional to the (5-function is usually neglected, 
and it does not contribute to the Purcell factor. How- 
ever, in hyperbolic media this singular term is complex 
even for vanishing losses, and it determines the diverging 
radiative decay rate. 
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Appendix B: Cartesian presentation of the field of 
dipole orthogonal to the anisotropy axis 



Cartesian components of Eq. ([26]) for the dipole ori- 
ented along y axis, perpendicular to the anisotropy axis 
z, read:[ 
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where 

r e = ^e\\(x 2 + y 2 ) + e±z 2 , 

r Q = \/e^(x 2 + y 2 + z 2 ). 

Magnetic field reads 

7 o . e lkre ( 1 \ y 2 z 

H x = -pk e\\y/e±—^- 1 



r 2 \ i/cr e / x 2 + y 2 
-pk 2 e ±v feI^ Y - 1 - — 

r 2 i ,,2A2 ' 



r r e J (x 2 + ?/ 2 ) 



n ikr e 



^IIV^Vl 1 -^)^ <B2) 
2 e*^ / 1 \ 



■ ipkyfs^ 



r \ ikr ) x 2 + ?/ 2 



r 

ikr Q 



(x 2 + ?/ 2 ) 2 ' 
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